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We determine the ability of future experiments to observe supersymmetric contributions to the rare decays 
B — > Xs'y and B — > X s l~*~l~ . A global fit to the Wilson coefficients which contribute to these decays is performed 
from Monte Carlo generated data. This fit is then compared to supersymmetric predictions for several different 
patterns of the superpartner spectrum. 



The first conclusive observation of penguin 
mediated processes, the exclusive B — > K*j and 
inclusive B — > X s j, by CLEO El has placed the 
study of rare B decays on new ground. These fla- 
vor changing neutral current (FCNC) transitions 
provide an essential opportunity to test the Stan- 
dard Model (SM) and offer a complementary strat- 
egy in the search for new physics by probing the 
indirect effects of new particles and interactions 
in higher-order processes. With the expected high 
luminosity of the B- Factories presently under con- 
struction (and the associated advanced detector 
technology), radiative B decays will no longer be 
rare events, and the exploration of FCNC transi- 
tions can continue by probing decay modes with 
even smaller predicted branching fractions. The 
cleanest rare decay which occurs at a rate acces- 
sible to these machines is B — » X s £ + £~ . In fact, 
experiments at e + e~ and hadron colliders are al- 
ready closing in on the observation □ of the ex- 
clusive modes B — * K^*H + £~ with £ = e and fx, 
respectively. Once this decay is observed, the uti- 
lization of the kinematic distributions of the £ + £~ 
pair, such as the lepton pair invariant mass dis- 
tribution and forward backwarcLasymmetrytl, and 
the tau polarization asymmetryQ in B — > X s t + t~ , 
together with B(B — > X s j) will provide a strin- 
gent test of the SM. In this talk we determine the 
ability of £?-Factories to probe possible supersym- 
metric contributions to these decays. 

Softly broken supersymmetry (SUSY) is a de- 
coupling theory, thus making it a challenge to 
search for its effects through indirect methods. 
However, a promising approach is to measure ob- 
servables where supersymmetry and the SM arise 
at the same order in perturbation theory. In this 
case the SUSY contributions do not suffer an ex- 
tra a/47r reduction compared to the SM ampli- 
tudes. The relative ratio between the lowest order 



SM amplitudes and those of supersymmetry could 
then be 0(1) if rh ~ Mw- Rare B-decays could 
then provide such an opportunity for discovering 
indirect effects of supersymmetry uu. 

The effective field theory for b — > s transitions 
which incorporates QCD corrections is governed 
by the Hamiltonian 
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where the Oi are a complete set of renormalized 
operators of dimension six or less which medi- 
ate b — > s transitions and are catalogued in, e.g., 
Ref. 7. The Ci represent the corresponding Wil- 
son coefficients which arc evaluated pcrturbatively 
at the electroweak scale where the matching con- 
ditions are imposed and then evolved down to the 
renormalization scale /i ~ mf,. The expressions for 
Ci(Mw) are given by the Inami-Lim functions^. 

For B — > X s £ + £~ this formalism leads to the 
physical decay amplitude (neglecting m s ) 
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where q 2 represents the momentum transferred to 
the lepton pair. We incorporate the NLO analysis 
for this decay which has been performed in Buras 
et al. u, where it is stressed that a scheme inde- 
pendent result can only be obtained by including 
the leading and next-to-leading logarithmic cor- 
rections to Cg(/i) while retaining only the leading 
logarithms in the remaining Wilson coefficients. 
The residual leading fi dependence in Cg(/x) is can- 
celled by that contained in the matrix element of 
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Table 1: Values of the Wilson coefficients for several choices 
of the renormalization scale. Here, we take = 4.87 GeV, 
m t = 1 75 GeV, and a s (M z ) = 0.118. 





jtt = TOft/2 


H = m b 


jJL = 2m b 




-0.371 


-0.312 


-0.278 


c 9 


4.52 


4.21 


3.81 


Cio 


-4.55 


-4.55 


-4.55 



09, yielding an effective value . The effec- 
tive value for Cy (fi) refers to the leading order 
scheme independent result. The operator Oio does 
not renormalizc. The numerical estimates (in the 
naive dimensional regularization (NDR) scheme) 
for these coefficients are displayed in Table [j]. The 
reduced scale dependence of the NLO versus the 
LO corrected coefficients is reflected in the devi- 
ations ACg(^) < ± 10% and AC 7 e// (/i) ~ ±20% 
as /i is varied. We find that the values of the co- 
efficients are much less sensitive to the remaining 
input parameters, with ACg^LACy^//) < 3%, 
varying a s (M*) = 0.118 ± 0.00300, and mf ys = 
175 ± 6 GeV O. The resulting inclusive branch- 
ing fractions (which are computed by scaling the 
width for B — > X s £ + £~ to that for B semi- 
leptonic decay) are found to be (6.25iJ'g 3 ) x 10 -6 , 
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e,/i, and r, respectively, taking into ac- 
count the above input parameter ranges, ms well 
as B s i = B(B -> Xiv) = (10.23 ± 0.39)%E3, and 
m c /m b = 0.29 ±0.02. 

The operator basis for the decay B — ■> X s ^f 
contains the first eight operators in the effective 
Hamiltonian of Eq. (0). The leading logarithmic 
QCD corrections have been completely resummed, 
but lead to a sizeable /i dependence of the branch- 
ing fraction and hence it is essential to include the 
next-to-leading order corrections. In this case, the 
calculation involves several steps, requiring NLO 
corrections to both Cy and the matrix element 
of 07. For the matrix element, this includes the 
QCD bremsstrahlung corrections Ej b — + S7 + g, 
and the NLO virtual corrections which have re- 
cently been completed Lj. Summing these con- 
tributions to the matrix elements and expanding 
them around /i = m b , one arrives at the decay 
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Here, the quantities 7^ are the entries of the ef- 
fective leading order anomalous dimension matrix, 
and the rt are computed in Greub et al. Eil, for 
i = 2, 7, 8. The first term in Eq. g, Cj ff {^), must 
be computed at NLO precision, while it is consis- 
tent to use the leading order values of the other 
coefficients. For Cy the NLO result entails the 
computation oiLthe 0(a s ) terms in the match- 
ing conditions t3, and the renormalization group 
evolution of Cy (/i) must be computed using the 
0{oq) anomalous dimension matrix. Preliminary 
NLO results for these anomalous dimensions have 
recently been reported 113, with the conclusion be- 
ing that in the NDR scheme the NLO correction 
to Cy (/i) is small. Therefore, a good approxima- 
tion for the inclusive width is obtained by employ- 
ing the leading order expression for Cy (/x) , with 
the understanding that this introduces a small in- 
herent uncertainty in the calculation. We then 
find the branching fraction (again, scaling to semi- 
leptonic decay) 

B(B -> X s j) = (3.25 ± 0.30 ± 0.40) x 10~ 4 , (5) 

where the first error corresponds to the combined 
uncertainty associated with the value of mt and 
fx, and the second error represents the uncer- 
tainty from a s (Mz), B s i, and m c /m b . This is well 
within the range observed by CLEO El which is 
B = (2.32±0.57±0.35) x 10" 4 with the 95% C.L. 
bounds of 1 x 10~ 4 < B(B -> X sl ) < 4.2 x 10~ 4 . 

Measurements of B(B — > X s j) alone con- 
strain the magnitude, but not the sign, of Cy (//). 
We can write the coefficients at the matching scale 
in the form d(M w ) = Cf M {M w ) + Cl iew (M w ), 
where Cf ew (Mw) represents the contributions 
from new interactions. Due to operator mixing, 
B — ► X s "f then limits the possible values for 
Cf ew {Mw) for i = 7, 8. These bounds are summa- 
rized in Fig. 0. Here, the solid bands correspond to 
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Figure 1: Bounds on the contributions from new physics to 
C73. The region allowed by the CLEO data corresponds 
to the area inside the solid diagonal bands. The dashed 
bands represent the constraints when the renormalization 
scale is set to fi = mi,. The diamond at the position (0,0) 
represents the standard model. 

the constraints obtained from the current CLEO 
measurement, taking into account the variation of 
the renormalization scale mj,/2 < /1 < 2rrib, as well 
as the allowed ranges of the other input parame- 
ters. The dashed bands represent the constraints 
when the scale is fixed to /i = mi,. We note that 
large values of C% ew (Mw) are allowed even in the 
region where C? ew (M w ) ~ 0. 

Measurement of the kinematic distributions 
associated with the final state lepton pair in B — > 
X s £ + £~ as well as the rate for B — > X s -f allows 
for the determination of the sign and magnitude 
of all the Wilson coefficients for the contributing 
operators in a model independent fashion. We 
have performed a Monte Carlo analysis in order 
to ascertain how much quantitative information 
will be obtainable at future B-factories and follow 
the procedure outlined in Ref. 5. For the pro- 
cess B — > X s £ + £~, we consider the lepton pair 
invariant mass distribution and forward-backward 
asymmetry for i = e, /a, t, and the tau polarization 
asymmetry for B — > X s t + t~. A three dimen- 
sional x 2 fit to the coefficients CV,9,iq(m) is per- 
formed for three values of the integrated luminos- 
ity, 3 x 10 7 , 10 s , and 5 x 10 s BB pairs, correspond- 
ing to the expected e + e _ £?-factory luminosities of 
one year at design, one year at an upgraded accel- 
erator, and the total accumulated luminosity at 
the end of the programs. The 95% C.L. allowed 
regions as projected onto the Cg(/i) — Cio(/i) and 
C7 (/x) — Cio(/i) planes are depicted in Figs. ||(a- 
b) , where the diamond represents the central value 
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(b) 95% C.L. 
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Figure 2: The 95% C.L. projections in the (a) Cg — Cio 
and (b) — Cio planes, where the allowed regions lie 

inside of the contours. The solid, dashed, and dotted con- 
tours correspond to 3 X 10 7 , 10 s , and 5 X 10 s BB pairs. 
The central value of the SM prediction is labeled by the 
diamond. 

for the expectations in the SM. We see that the 
determinations are relatively poor for 3 x 10 7 BB 
pairs and that higher statistics are required in or- 
der to focus on regions centered around the SM. 

We analyze the supersymmetric contributions 
to the Wilson coefficients in terms of the quantities 

= cr y (M W ) _ = cr w (M W ) 

Ul ~ Cf M (M w ) - Cf M (M w ) ' W 

where C^ usy (Mw) includes the full SM plus su- 
perpartner contributions. Ri is meant to indicate 
the relative fraction difference from the SM value. 

Supersymmetry has many potential sources 
for flavor violation. The flavor mixing angles 
among the squarks are a priori separate from the 
CKM angles of the SM quarks. We adopt the 
viewpoint here that flavor-blind (diagonal) soft 
termsE] at the high scale are the phenomenologi- 
cal source for the soft scalar masses, and that the 
CKM angles are the only relevant flavor violat- 
ing sources. The spectroscopy of the supersym- 
metric states is quite model dependent and here 
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we analyze two possibilities. The first is the fa- 
miliar minimal supergravity model. In this case 
all the supersymmetric states follow from a com- 
mon scalar mass and a common gaugino mass at 
the high scale. The second possibility is to relax 
the condition of common scalar masses and allow 
them to take on uncorrelated values at the low 
scale while still preserving gauge invariance. 

We begin by searching over the full parame- 
ter space_of the minimal supergravity model. We 
generateliS these models by applying common soft 
scalar and common gaugino masses at the bound- 
ary scale. The tri-scalar A terms are also input at 
the high scale and are universal. The radiative 
electroweak symmetry breaking conditions yield 
the B and /x 2 terms as output, with a sign(/x) 
ambiguity left over. (Here /i refers to the Hig- 
gsino mixing parameter.) We also choose tan /3 
and restrict it to a range which will yield pertur- 
bative Yukawa couplings up to the GUT scale. We 
have generated thousands of solutions according 
to the above procedure. The ranges of our input 
parameters are < m < 500 GeV, 50 < my 2 < 
250 GeV, -3 < A /m < 3, 2 < tan/3 < 50, and 
we have taken m t = 175 GeV. Each supersym- 
metric solution is kept only if it is not in viola- 
tion with present constraints from SLC/LEP and 
Tevatron direct sparticle production limits, and it 
is out of reach of LEP II. For each of these re- 
maining solutions we now calculate Rt-iq. Our 
results are shown in the scatter plots of Fig. |^ in 
the (a) Rj — R$ and (b) Rg — R w planes. The 
diagonal bands represent the bounds on the Wil- 
son coefficients as previously determined from our 
global fit. We see that the current CLEO data on 
B — > X s j already places signigicant restrictions 
on the supersymmetric parameter space. 

The first thing to note from the figure is that 
large values of i?7 and R$ are generated, and that 
they are very strongly correlated. These large ef- 
fects arise from models with \fi\ < 400 GeV. This 
is because light charged Higgsinos (or rather light 
charginos with a large Higgsino fraction) are nec- 
essary in order to obtain a large effect on the Wil- 
son coefficients. We see that the values of Rg and 
Rio are bounded by about 0.04, which is small 
compared to the range for R7. The main rea- 
son for these smaller values is the dependence on 
the bottom Yukawa Af, oc l/cos/3. i?7 also has 
a contribution directly dependent on this 1 / cos /3 
Yukawa factor, however the other multiplicative 



: 1 






1 l' 
1 1 


1 1 1 1 ; 1 " : 








1 1 
■■ii. - 1 1 


/ \ 

/ 


1 1 






' 1 1 

1 1 
1 1 
1 \ 
\ \ 

, ih 





-7.5 -5-0 -2-5 0.0 2.5 5.0 



R 7 



- ' ' ,'. ' 1 ' 1 ' 1 1 1 ' 1 ' 1 ' 1 ' 1 - 
SJvrS^^i — 


US 












"X' 








■■■'.•!}>■''■,■ 





0.00 0.01 0.02 0.03 0.04 



R 9 

Figure 3: (a) Parameter space scatter plot of R7 vs. ijg in 
the minimal supergravity model. The allowed region from 
CLEO data lies inside the solid diagonal bands. The dashed 
band represents the potential 10% measurement from the 
previously described global fit to the coefficients, (b) Pa- 
rameter space scatter plot of Rg vs. -Rio- The global fit to 
the coefficients obtained with 5 X 10 s BB pairs corresponds 
to the region inside the diagonal bands. 

terms associated with Af, are the large top Yukawa 
and a large kinematic loop factor. Rg and i?io 
do not have such additional factors due the chiral- 
ity structure of these operators and the require- 
ment that leptons and sleptons only couple via 
SU(2) and (7(1) gauge couplings. These con- 
ditions, along with the correlations between the 
mass spectra dictated by minimal supergravity re- 
lations, render the minimal supergravity contribu- 
tions to i?9.io essentially unobservable. 

We now adopt our second, more phenomeno- 
logical, approach. The maximal effects for the pa- 
rameters Ri can be estimated for a superparticle 
spectrum independent of the high scale assump- 
tions. However, we still maintain the assump- 
tion that CKM angles alone constitute the sole 
source of flavor violations in the full supersym- 
metric lagrangian. We will focus on the region 
tan/3 < 30. The most important features which 
result in large contributions are a light t\ state 
present in the SUSY spectrum and at least one 
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Figure 4: The maximum value of (a) _Rg and (b) Rio 
achievable for general supersymmetric models. The top 
solid line comes from t — contribution and is displayed 
versus the mass. The bottom solid line is from ti — x^f 
contribution with tan f) = 1 and is shown versus the Xi 
mass. The dashed line is the ti — xf 1 contribution with 
tan /3 = 2. The other mass parameters which are not plot- 
ted are chosen to be just above LEPII and Tevatron's reach. 



light chargino state. For the dipole moment oper- 
ators a light Higgsino is most important. A pure 
higgsino and/or pure gaugino state have less of an 
effect than two mixed states when searching for 
maximal effects in Cg and Cio and we have found 
that Mi ~ 2/i is optimal. 

Fig. H displays the maximum contribution to 
i?9 i io versus an applicable SUSY mass scale. The 
other masses which are not plotted (ti, Z^, etc.) 
are chosen to be just above the reach of LEPII or 
the Tevatron, whichever gives better bounds. We 
see that the maximum size of -R9.10 is much larger 
than what was allowed in the minimal supergrav- 
ity model. This is due to the lifted restriction 
on mass correlations. Light sleptons, sneutrinos, 
charginos, and stops are allowed simultaneously 
with mixing angles giving the maximal contribu- 
tion to the i?i's. However, we find that the maxi- 
mum allowed values for -R^io are still much less 
than unity. Earlier we determined that B fac- 



tory data would be sensitive to ARg > 0.3 and 
Ai?io > 0.08 at the highest luminosities, and so 
the largest SUSY effect would give a l — 2a signal 
in i?g,io, hardly enough to be a compelling indica- 
tion of physics beyond the standard model. 

Given the sensitivity of all the observables it 
is instructive to narrow the focus to Cr(Mw)- 
There exists the possibility that one eigenvalue 
of the stop squark mass matrix might be much 
lighter than the other squarks due to the large top 
Yukawa and the mixing term mt (A t — [i cot /3) in 
the stop mass matrix. We then present results for 
Ct(Mw) in the limit of one light squark, namely 
the ti, and light charginos. We allow the t\ to 
have arbitrary components of and tu since cross 
terms can become very important. This is es- 
pecially noteworthy in the high tan/3 limit. We 
note that the total supersymmetric contribution 
to Ci(M\\r) will depend on several combinations 
of mixing angles in both the stop and chargino 
mixing matrices and cancellations can occur for 
different signs of fitB. 

The first case we examine is that where the 
lightest chargino is a pure Higgsino and the light- 
est stop is purely right-handed: \i ~ H ± , t\ ~ 
ti?.. The resulting contribution to is shown as 
a function of the Ir mass in Fig. || (dashed line) 
for the case of chargino masses out of reach of 
LEP II (m + > M w ). Note that the SUSY con- 
tribution to Ct(Mw) in this limit always adds con- 
structively to that of the SM. Next we examine the 
limit where the only light chargino is a pure Wino. 
The effects of a light pure Wino are generally small 
since (i) it couples with gauge strength rather than 
the top Yukawa, and (ii) generally supersymmetric 
models do not yield a light tt necessary to couple 
with the Wino. This contribution to R-? is shown 
in Fig. H (dotted line). As expected, we see that 
this contribution is indeed small. Now we discuss 
our third limiting case. As mentioned above, what 
we mostly expect in minimal supergravity models 
is a highly mixed t\ state; here, for large effects, 
it is crucial that there be substantial tn and 
contributions to t\. We find that in this case large 
tan /3 solutions (tan /3 > 40) can yield greater than 
0(1) contributions to R7 even for SUSY scales of 
1 TeV! Low values of tan (3 can also exhibit sig- 
nificant enhancements. This is demonstrated for 
tan/3 = 2 in Fig. || (solid line). Note that in this 
case large contributions are possible in both the 
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Figure 5: Contributions to Ft? in the different limits de- 
scribed in the text. The top solid line is the charged H^/t 
contribution versus m H ± . The bottom solid line is the 
Xi ft\ contribution versus TTis± where both the chargino 
and stop are maximally mixed states with /i < 0. The 
dashed line is the H^/Ir contribution, and the dotted line 
represents the W^/t\ contribution. These two lines are 
both shown as a function of \^ mass. All lines are for 
tan (3 = 2 and rat = 175 GeV. We have set all other masses 
to be just above the reach of LEPII. 



negative and positive directions of R-j depending 
on the sign of fi. We note that this is a region of 
SUSY parameter space which is highly motivated 
by 50(10) grand unified theories. 

In this talk we have studied the effects of su- 
persymmetry to the FCNC observables concern- 
ing b — ► s transitions, and we have seen that de- 
viations from the standard model could be de- 
tected with supersymmetric masses even at the 
TeV scale. This is especially true if tan/3 is very 
high. The large tan/3 enhancements in the b — > s 
processes are unique; it is thus possible that the 
first distinct signs of supersymmetry could come 
from deviations in rare B decays. One, of course, 
would like direct confirmation of such a deviation, 
if observed at i?-Factories, and collider programs 
could provide it. 
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